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The near edge structure (XANES) in K-edge X-ray absorption spectroscopy (XAS) is a widely 
used tool for studying electronic and local structure in materials. The precise interpretation of 
these spectra with the help of calculations is hence of prime importance, especially for the study 
of correlated materials which have a complicated electronic structure per se. The single particle 
approach, for example, has generally limited itself to the dominant dipolar cross-section. It has 
long been known however that effects beyond this approach should be taken into account, both due 
to the inadequacy of such calculations when compared to experiment and the presence of shake-up 
many-body satellites in core-level photoemission spectra of correlated materials. This effect should 
manifest itself in XANES spectra and the question is firstly how to account for it theoretically 
and secondly how to verify it experimentally. By using state-of-the-art first principles electronic 
structure calculations and Is photoemission measurements we demonstrate that shake-up many- 
body effects are present in K-edge XAS dipolar spectra of NiO, CoO and CuO at all energy scales. 
We show that shake-up effects can be included in K-edge XAS spectra in a simple way by convoluting 
the single-particle first-principles calculations including core-hole effects with the Is photoemission 
spectra. We thus describe all features appearing in the XAS dipolar cross-section of NiO and CoO 
and obtain a dramatic improvement with respect to the single-particle calculation in CuO. These 
materials being prototype correlated magnetic oxides, our work points to the presence of shake-up 
effects in K-edge XANES of most correlated transition metal compounds and shows how to account 
for them, paving the way to a precise understanding of their electronic structure. 



I. INTRODUCTION 

Low-energy excitations of correlated materials arc 
strongly influenced by electron-electron interaction ef- 
fects beyond the single particle approximation. This is 
particularly evident in core-level photoemission spectra 
(XPS) of transition metal compounds where the occur- 
rence of many-body satellites is well documented (for a 
review see Ref. [D). As for x-ray absorption, interpre- 
tation of the dipolar K-edge XAS cross-section heavily 
relies on standard single-particle first principles calcula- 
tions— _ — that neglect shake-up excitations. Since dipolar 
1,2,3 XAS mostly samples d-states of the absorbing atom 
which are more prone to effects of correlation than p- 
states, one normally assumes that shake-up effects are 
visible mostly in 1^3 XAS and not in K XAS. However a 
recent work& shows that in NiO the single particle dipo- 
lar K-edge spectrum misses some near-edge and far-edge 
features present in the experimental measured one. 

We concentrate on shake-up many body excitations 
arising from a valence electron excitation following the 
creation of a core hole by the incident x-ray 7 -. In the past, 



the description of shake-up effects in core-hole photoe- 
mission spectra has been investigated in the framework 
of quantum-chemical calculations—, by using approaches 
based on model hamiltoniansi^ii^ or first-principles mod- 
ified approaches 11 . The occurrence of these effects in 
XPS is well established but they have also been shown to 
occur in M4.5 edges of mixed-valent compounds^ and in 
L/2,3 X-ray absorption spectra (XAS) of transition met- 
als and rare earths compound o 13 ' 14 and were proposed 
as a possible explanation of the double peak structure in 
dipolar K-edge XAS of high T c cuprates^ and copper 
compounds in general^. However this attribution was 
questioned in Ref^i^ 7 - and the double peak structure was 
suggested to be single particle in origin. 

Nailing down the importance of these effects has been 
difficult due to complications related to many-body cal- 
culations but also to the paucity of experimental Is pho- 
toemission spectra. In this work, following earlier sugges- 
tions, we demonstrate that shake-up manybody effects 
can be included in a simple way in K-edge XAS spectra 
by convoluting the single particle first principles calcula- 
tions with experimental Is photoemission spectra, some 
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of which we have freshly measured. We show that this 
procedure explains all features in K-edge XAS spectra 
of NiO and CoO and strongly improves the agreement 
with experimental data in CuO. Our work points out the 
relevance of these effects in K-edge dipolar XAS of all 
compounds displaying multiple structures in photoemis- 
sion spectra. 

The structure of the paper is the following. In section 
III A[ following Ref H3, we briefly sketch the demonstra- 
tion of the convolution formula relating the many body 
XAS cross section to the single particle one via the pho- 
toemission spectrum. We then present experimental de- 
tails concerning our measured photoemission spectra in 
sec. IIII1 The experimental results and the theoretical 
undestanding are presented in sec. IIVI 



II. THEORY 

A. Shake-up theory 

Shake up satellites are many-body peaks present in 
core-electron spectra. They originate from a valence elec- 
tron excitation following the creation of a core hole by the 
incident x-ray-£. Quantum chemical calculations of shake- 
up satellites have been recently reviewed by Carravetta 
and Agren^. 

An electric dipole transition between two Slater de- 
terminants built from the same set of orbitals do not 
allow for shake-up satellites. Indeed, the orthogonality 
of orbitals allows for only one transition from the core 
level to the empty one. Therefore, a shake-up can only 
be obtained by describing the (initial) state with a lin- 
ear combination of Slater determinants or by using dif- 
ferent orbitals for the initial and final determinants^. 
The first approach was extensively used by Sawatzky and 
collaborators^ . Here we use dipole transitions between 
single Slater determinants using non-orthogonal orbitals, 
the orbitals of the final state being relaxed in the presence 
of the core hole. 

The possibility of describing the electronic state of NiO 
by a single Slater determinant was suggested in Refs 
1 9~ 21 . Moreover, relaxed Slater determinants can some- 
times describe a state much better than the sum of a 
small number of unrelaxed Slater determinants 1 ^. 

A single Slater determinant is also the non-interacting 
ground state of the Kohn-Sham version of density fun- 
tional theory (DFT). The corresponding Kohn-Sham or- 
bitals are usually considered to have no physical mean- 
ing. This would be a problem for our approach that 
calculates electric dipole transitions between these or- 
bitals. The success of DFT calculations of XAS seems to 
indicate that Kohn-Sham orbitals are physically mean- 
ingful and indeed Gidopoulos 22 discovered that the non- 
interacting Kohn-Sham ground state is the best approx- 
imation of the true ground state in a subtle way. To de- 
scribe his finding, let h(r) = — h 2 A/2m + v(r) be a one- 
body potential and H v = J2i M r ») be the correspond- 



ing non-interacting many-body Hamiltonian. Denote by 
\iffv) the (Slater determinant) ground state of H v and 
by |^) the ground state of the interacting Hamiltonian 
H. By the Rayleigh-Ritz minimum principle we have 
(*|jff u |*) - (*„|flo|* ) > 0. Gidopoulos proved that 
the potential v that minimizes this difference is precisely 
the Kohn-Sham potential. In that sense, the Kohn-Sham 
determinant and the Kohn-Sham potential provide the 
best single-particle description of the ground state of an 
interacting system. 

Therefore, it is relevant to describe shake-up processes 
with non-orthogonal Slater determinants. Other calcula- 
tions were carried out within this framework by Tyson 2 ^, 
who could calculate double-electron excitations in XAS 24 
for LN4.5-edges. Similar calculations for x-ray photoe- 
mission spectroscopy are more common 1 ^. 



B. Cross section 

The manybody X-ray photoemission cross section can 
be written as 2 ^ 

o_ 

axP S (e k ) = T J2\(k^f(N~l)\T\MN))\ 2 
f 

6(e k -fuv-Ei{N) + E f (N-l)) (1) 

where e k is the photoelectron kinetic energy, Ei (N) is the 
energy of the N electrons ground state |$;(iV)), Ef(N — 
1) and \^j(N — 1)) characterize the excited energy and 
state of the N—l electron system with a core hole and haj 
is the energy of the incident X-ray beam. The electric- 
dipole transition operator is labeled T. The transform 
Ixps(t) of the XPS cross-section is defined as, 

p+oo 

a X p S (e) =2Re dte u +*I X ps(t) (2) 
Jo 

where e + = e + ii~i and Eq. [2] has to be understood as the 
limit for rj — > + . 

The manybody X-ray absorption cross section in the 
dipolar approximation can be written as 

vxasH = ^£|<*/(A0|W1$ 4 (A0)| 2 
/ 

S(E f (N) - Ei(N) - tiu) (3) 

where now uj is the energy of the incident beam and W is 
proportional to the dipole transition operator M, namely 
W = ^/2irh 2 ujaoM. 

Similarly to what was done for the case of XPS and 
using a similar notation, we can define the the transform 
IxAs{t) of the XAS cross-section as 

(T X As(t) =2 Re/ dte l ^ l I X As{t) (4) 
Jo 
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Under the assumption that both &i and tyj ps are sin- 
gle determinant states, Ohtaka and Tanak o 26 i demon- 
strated that 

I X AS(t) = Ixps{t)h(t) (5) 

where both the terms Ixps{t) and Io(t) (see Eq. 4.43 
in Ref. 27 ) includes manybody shake-up processes at all 
orders. Eq. [5] holds for a generic static core-hole poten- 
tial. A similar relation was found for the less-general 
case of a contact core-hole potential by Nozieres and 
DeDominicis 28 using the linked cluster theorem. 

If shake-up processes arc neglected only in the Io(t) 
term then I (t) reduces to IxAs(t)i name ly the trans- 
form of the single-particle XAS cross section calculated 
in the presence of a static core-hole potential. Thus, it is 
possible to include to some extent many body effects in 
the XAS cross section o~xas{u) by performing the convo- 
lution between the measured XPS cross section crxps( e ) 
that fully includes manybody shake-up processes and 
the single-particle calculated XAS cross section o- S x AS (ui), 
namely 

ctxas(u) = J d e ^xps( e ) u xAs( u - £ ) ( 6 ) 

In our work the single-particle cross section o~x AS (uj) 
is calculated in the framework of density functional the- 
ory with inclusion of a static core- hole and (J^p S {w) is 
measured. 



C. Technical details 

The single particle XAS cross section (T^ AS {(J) is cal- 
culated in the framework of density functional theory 
using the implementation of Ref. 2 distributed with the 
Quantum-Espresso^ distribution. The technical details 
for the NiO calculation are the same as in ref. 6 . We 
used norm-conserving pseudopotentials with inclusion of 
semicore states. The energy cutoffs used in the calcu- 
lations were 140 Ryd and 160 Ryd for CuO and CoO, 
respectively. In the case of CoO, we neglect the tetrag- 
onal structural distortion below the 290K Neel temper- 
ature and adopt magnetic and crystal structures similar 
to those of NiO. The electron- momentum grids for the 
Brillouin zone integration and the choice of the supercell 
for the XAS calculation are the same as for the NiO case 
in ref£. The CuO XAS cross-section was calculated in 
the supercell obtained by doubling the antiferromegnetic 
cell along the shortest direction. The antiferromagnetic 
cell is obtained from the non-magnetic one by defining as 
new lattice vectors a' = a + c, b' = b, and c' = a c 
where a,b, and c are the direct lattice vectors. We then 
use a 3 x 3 x 3 electron-momentum grid in the supercell 
to obtain the self consistent charge density and a 3 x 3 x 3 
electron-momentum grid in the supercell to calculate the 
XAS cross-section. 
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FIG. 1: figure 1: Experimental (circle) and fitted (lines) Is 
photoemission spectra in NiO 

Finally we employ the DFT+U approximation in all 
case with U — 7.75 eV and U = 11.1 eV for CoO and 
CuO respectively. These values of the Hubbard repulsion 
are calculated from first principles using the method of 
ref. 30 . 



III. EXPERIMENT 

The Ni-ls photoemission spectrum in NiO were mea- 
sured at the HIKE station of the KMC-1 beamline at 
BESSY- 3 ^. The spectra were recorded with a SCI- 
ENTA R4000 photoelectron analyzer placed at 90° from 
the x-ray beam. The incident x-ray beam (9 keV) was 
monochromatized by a pair of Si(422) crystals providing 
^500 meV energy bandwidth. To avoid charging effects, 
a 25 nm thick NiO thin film was grown on a Ag sub- 
strate in the presence of oxygen, and capped by 3 nm of 
MgO. The growth of NiO was found fully epitaxial with 
the NiO(OOl) direction parallel to Ag(001) as confirmed 
by the LEED patterns. The sample was positioned at a 
grazing angle of 89.99° from the incident x-rays in order 
to reduce the penetration depth of photons and enhance 
the photoelectron yield. 

The XAS spectra of NiO and CoO were borrowed from 
Refs^ and^i respectively. 

IV. RESULTS 

A. Nickel Oxide 

The measured Is photoemission spectra of NiO are 
shown in Fig. [T] The fit to the data is consistent with 
a three peak structure in the 590-610 eV energy region. 
The results closely resemble 2p 3 / 2 Ni photoemission in 
NiO 10 -^ in this energy region. In literature the attri- 
bution of the different features in 2p 3 / 2 Ni XPS is very 
controversial and was subject to several reinterpretations. 
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FIG. 2: Convolution of experimental XPS (this work) and 
single particle XAS calculation in NiO at the Ni K-edge. Con- 
volution with two or three component is shown with solid line 
and dashed line respectively (middle spectra). Experimental 
XAS data are from Refi^. 



Van Veenendaal and Sawatzky 9 attributed the main fea- 
ture at high energy to a 2p 5 3d 9 L state, where L means a 
hole in the ligand state. The satellite of the main peak 
(shoulder) was attributed to non-local screening coming 
from the nearest neighbours Ni atoms, while the lower en- 
ergy satellite at « 596 eV was attributed to a 2p 6 3d 10 L. 
Recently this was reconsidered in ref.— where the main 
feature was attributed to a 2p 5 3d 9 Z, where Z is a Zhang- 
Rice k-dispersing bound stated. The shoulder of the 
main peak is attributed to a 2p 5 3d 9 L state and the low- 
est energy feature to a 2p 5 3d 8 state. Here we show that, 
regardless of their attribution, the features measured in 
Is Ni NiO XPS are present also in the dipolar Ni K-edge 
XAS spectrum of NiO. 

In Fig. [2] we show the measured and calculated XAS 
cross sections. The single particle cross section is gen- 
erally in good agreement with the dipolar part of the 
measured spectrum except for the two peaks indicated 
by the letters F and H that are missing in the single par- 
ticle calculation. In order to determine if the missing 
excitations are manybody in nature, and eventually due 
to multi-determinant or shake-up processes, we then pro- 
ceed by using Eq. [B] and obtain new XAS spectra. We 
first perform the convolution using the complete three- 
peaks structure of the XPS spectra. We find that the con- 
volution of the DFT calculated XAS cross-section with 
the photoemission spectra greatly improves the agree- 
ment with experiments. In particular the missing peaks 
are now present in the spectrum and a better agreement 
occurs at all energy scales. The F and H peaks are then 
replicas of the single particle E and G peaks respectively 
and are manybody in nature. We can further test to 
what extent this interpretation is robust by altering the 
XPS spectrum before convolution with the theoretical 
single-particle XAS calculation. We do this for NiO by 
artificially supressing the main peak of the XPS spec- 



trum, leaving the shoulder and the satellite. We find 
that the resulting absorption spectrum (dashed line Fig 
xxx)is less in agreement with the experimental spectrum 
supporting our interpretation. 



B. Cobalt Oxide 
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FIG. 3: Convolution of experimental XPS data of Refs.— ^ 
and single particle XAS calculation in CoO at the Co K-edge. 
Experimental XAS data are from Ref.—. 

Co Is photoemission data of CoO are not available in 
literature. However 2p 3 / 2 and 3s Co XPS dat a 35 ' 37 are 
extremely similar and composed by two main peaks and 
a small shoulder at low energy visible only in the 3s data. 
We then consider 3s photoemission data of Ref. 3 - and 
fit them with a two-peak structure and neglect the very 
small shoulder, invisible in 2p photoemission data. The 
results are shown in Fig. [3] 

The situation is very similar to NiO, namely the F peak 
is missing missing from the single-particle spectrum and 
the H peak is weak. Convolution with photoemission 
improves substantially the agreement although the main 
edge peak is narrower than the experimental data. 



C. Copper Oxide 

The copper oxide CuO has a monoclinic crystal struc- 
ture with symmetry group C/2c . The dipolar part in 
CuO was measured in refi 38 ' 39 . We follow the notation 
of refs i 38 ' 39 and label the configuration of the crystal with 
respect to the incident beam by the three angles (8, 0, ip). 
These angles correspond to the three rotation angles of 
the goniometer. In particular, when the three angles are 
zero, then the polarization is parallel to the 9 axis of the 
goniometer and to the c-axis of the crystal. At zero an- 
gles, the plate holding the sample is orthogonal to the 
incident beam and parallel to one of the plaquette chains 
in the crystal (see Reb^ 9 - Fig. 10 for more details). 
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FIG. 4: Convolution between measured XPS from Ref.— and 
single particle XAS calculation at the Cu K-edge. Experi- 
mental XAS data are from 3 ^ 



Given the low symmetry of the crystal, the polariza- 
tion dependence of CuO K-edge XAS spectra is very 
complicated, as can be seen in Fig. |4j The calculated 
single particle spectra are in strong disagreement with 
experiments. Both the peak positions and the polar- 
ization dependence of the intensities disagree with the 
measured data. In order to see if the disagreement is 
due to the lack of manybody effects in the XAS cross- 
section, we consider the convolution with Cu Is photoc- 
mission spectral. Cu Is photoemission spectra of CuO 
are composed of two peaks, usually attributed to 3d 9 
and 3d 10 L. Performing the convolution with the calcu- 
lated single particle Cu K-edge XAS leads to an impres- 



sive improvement. The convoluted spectrum is in much 
better agreement with experimental data, demonstrating 
that the Cu K-edge XAS spectrum in CuO is dominated 
by shake-up manybody processes. 



V. CONCLUSIONS 

We have demonstrated that shake-up processes occur 
in dipolar K-edge XAS spectra of NiO, CoO, CuO. As 
these are prototype correlated transition metal oxides, 
we expect these excitations to be present in all XAS data 
of correlated materials. To be more precise, whenever 
charge transfer satellites occur in XPS core-hole spec- 
tra, then shake-up satellites must also occur in the cor- 
responding X-ray absorption edge, at all energy scales. 

We have also shown that a practical way to include 
these effects in first principle calculations is to perform 
the convolution with the XPS spectrum at the same edge, 
as suggested by Eq. [SI Despite the fact that Eq. [5] was 
obtained many years ago2£r— , we are currently unaware 
of other works trying to explicitly apply this equation to 
K-edge XAS by using state of the art calculations. Our 
work that full includes core-hole attraction and Hubbard 
U at the DFT+U level demonstrates that this approach 
is feasible and allows, for the first time, the attribution 
of all dipolar peaks in NiO and CoO dipolar K-edge XAS 
spectra. 

In Eq. [6J we neglected the additional many-body 
term o 26 ' 27 that are present in Io(i). These terms seem 
to be neglible in NiO and CoO, but could explain the re- 
maining discrepancy between theory and experiment in 
CuO. Further work is required to calculate these many- 
body correction terms. 
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